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Abstract 

We present a general discussion of the mean field dynamics of 
finite nuclei prepared under extreme conditions of temperature and 
pressure We compare the prediction of semi-classical approximation 
with complete quantum simulations. Many features of the dynamics 
are carefully studied such as the collective expansion, the evaporation 
process, the different time-scale... This study points out many quan- 
titative differences between quantum and semi-classical approaches. 
Part of the differences are related to numerical features inherent in 
semi-classical simulations but most of them are a direct consequence 
of the non treatment of nuclei as quantal objects. In particular, we 
show that because of a too strong damping in semi-classical approaches 
the expansion of hot nuclei is quenched and the speed of the collective 
motion reduced. 
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1 Introduction 



Semi-classical methods are widely used to describe the nuclear dynamics. In 
particular, considerable progresses have been made in semi-classical transport 
theory of Heavy Ion Collisions with the simulation of Boltzmann type, in 
which an averaged Pauli-blocked collision term is added on top of a mean- 
field evolution^ P], 0, H|. More recently, following the Boltzmann-Langevin 
(BL) theory, a stochastic term has been also introduced to take into account 
the strong fluctuations of observables ||. These approaches (BUU and BL) 
are of great interest not only because they are numerically tractable but also 
because they seem to be globally in good agreement with experimental data. 

As far as quantum features are concerned, semi-classical approaches es- 
sentially take into account the exclusion principle through the Pauli blocking 
in nucleon-nucleon collisions. Thus, many properties of the nucleus which 
are due to the quantum nature of the nuclear systems |10||, might not be 
correctly described in semi-classical approximations. For example, dereal- 
ization of wave-packets, wave dynamics, barrier transmission... are absent 
from semi-classical approaches, it is thus of great interest to check if the loss 
of these quantal features of nucleons does not change the evolution of a highly 
excited nucleus. The validity of semi-classical theories was already discussed 
several times and some qualitative differences have already been pointed- 



out |TTJ. However, we want to present a systematic quantitative comparison 
for the desexcitation of nuclei under extreme conditions. 

In this paper, we compare the evolution of a hot and compressed (or di- 
luted) nucleus in quantal calculation with the results of a semi-classical sim- 
ulation. Our aim is to bring quantitative information about the possible role 
of the quantal nature of the nucleons and about the accuracy of semi-classical 
approaches. Since, most of the discussed phenomena are already present at 
the mean-field level, we will restrict ourselves to a comparison between time 
dependent Hartree-Fock (TDHF) methods and its semi-classical counterpart: 
the Vlasov dynamics. Moreover, we will consider a simple scenario of the evo- 
lution of a hot spherical compound nucleus and we will analyze evaporation 
of particles, dynamics of collective motion, time-scales... 

In section 2, we present the considered scenario together with the various 



1 These approaches are often referred under different acronyms BUU, VUU, BNV[|[ [| 
0, ||... In this article we will use BUU (for Boltzmann-Ueling-Ulenbeck) which is the one 
initially proposed by Bertsch et al [Q, |[ |) 
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ingredients entering in the mean-field approximation. We will see that some 
of the differences between quantum and semi-classical simulations could be 
related to the finite sampling of the phase-space density, necessary in semi- 
classical methods. This numerical sampling of the phase-space generates a 
spurious finite range term in the semi-classical mean-field which was not in 
the original quantal description. This term modifies the properties of a finite 
piece of nuclear matter. When the quantum potential has a finite range 
component one may try to reduce the spurious effects by artificially reducing 
the range of the interaction used in the semi-classical simulations. We will 
show that this method is hardly a good approximation especially due to 
the density dependence of the various term in the interaction. In order to 
disentangle the role of this numerical artifact from a real difference between 
a quantal and a semi-classical treatment, a third type of simulation will be 
displayed: we will consider a TDHF simulation where a surface term has 
been added to the interaction in order to mimic the semi-classical numerical 
range. 

In section 3, we present results concerning the expansion dynamics of 
initial excited nuclei with different compression or dilatation. 

As expected, the semi-classical simulation seems to be a better approxi- 
mation of the TDHF simulation where the surface term has been introduced. 
Nevertheless, many differences remain, showing limits of semi-classical calcu- 
lations. For instance, evaporation rate, dynamics of the collective motion and 
associated time-scales are different. The condition under which a quantum 
system will expand toward low density regions will be finally discussed. 

2 A scenario to compare mean-field approaches. 

A good benchmark for a comparison between semi-classical and quantum dy- 
namics is provided by the model case introduced by D. Vautherin et al|13|. 
Indeed, they consider an excited spherical nucleus at various initial tempera- 
ture and with various density profiles. Then, they follow the evolution in the 
TDHF framework. The argument in favor of the use of a mean-field evolution 
is that the nucleon-nucleon collisions are not important for the dynamics of 
such thermalized system. Indeed, in the case of an equilibrated expanding 
source the gain and the loss terms in the collision part of the BUU approach 
almost cancel. 
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2.1 The interaction 



We will study the self-consistent mean field evolution of a spherical nucleus 
under extreme conditions. In the first part of this article we shall consider a 
density dependent effective two-body interaction: 

V(l, 2) = t 6& - # 2 ) + 1 3 p ( rj ^A 6{A - f 2 ), (1) 



which is a simplified version of the more general Skyrme force, but has been 
shown to be sufficient for the study of the expansion dynamics ||13||. In the 



mean-field approach, for spin-isospin saturated systems, this interaction cor- 
responds to a local potential energy 

E[p}(r,t) = ^ oP 2 (r,t) + ^P 2+<T M), (2) 

Note that, in the second part of this article, we shall add more terms to 
this interaction in order to take into account the finite range of the nuclear 
forces. Many of the useful results and derivations considering mean-field for 
a Skyrme force are derived in ref. ||14||. 

For our numerical applications, we have taken t = — lOOOMeVfm 3 , 
£3 = 15000MeVfm 6 and a = 1 (see ref.|13[). In the infinite nuclear mat- 
ter, this potential is globally in agreement with some more complicated 
parametrization ||15|| : this particular force yields to the following saturation 



point properties: E/A = —14MeV, a saturation density p = 0.1 5/m~ 3 and 
a compression modulus = 336.2MeV . This value corresponds to a rather 
stiff equation of state. Another interesting quantity is the "zero temperature 
spinodal density" , this is the density where the compressibility drops to 
zero at zero temperature. For the considered force, this density is equal to 
O.lfm' 3 . 



2.2 Constrained mean-field initialization 

In simulation of mean field type, the usual way for preparing systems at 
different temperatures and with different shapes -< Q y, where Q is the 
observable measuring the deformation, is to perform a constrained mean- 
field calculation with a constraining field XQ. For instance, Q can be the 
square radius operator and A is computed in order to obtain the desired 
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square radius -< f 2 >- or equivalently the desired compression factor r\ = 
(-< r 2 >~ I -< f 2 y) 3 ^ 2 (where -< f 2 ^ is the RMS of the ground state at 
zero temperature). We will also define the compression using the average 
density (p) in a sphere of radius 2fm around the center of the nucleus. The 
value of 2 fm is small enough compared to the RMS to avoid the bias of the 
surface and big enough to take into account not only the central s-orbitals 
but also others in order to focus on collective effects. 

The temperature can be either fixed arbitrary to study isothermal initial 
conditions or can be defined in order to get a given energy -< H y. This 
method is well adapted to generate isothermal initial conditions however 
when the nucleus is strongly compressed, the self-consistency of the mean- 
field may lead to some shape modification such as the apparition of hollow 
structures. Moreover, it appears rather difficult to generate diluted nuclei 
since negative A values leads to unbound systems. 

2.3 Isoentropic scaling initialization 

In order to overcome the limitations of the constrained mean-field approach 
and to be able to reach both high or low densities we have considered an 
alternative procedure based on the scaling of the nucleus. Since the scaling 
preserves the entropy, it is more natural to consider isoentropic families of 
initial nuclei. 

In order to define the initial density before scaling, we can use the con- 
strained mean field method described above with a small constraint param- 
eter A introduced to correctly treat the continuum (A = 0.25MeV fm~ 2 in 
our simulation). The initial temperatures T are defined in order to get the 
desired entropies S. 

For numerical reasons, we have chosen to neglect the temperature depen- 
dence of the mean-field potential in the solution of the constrained mean-field. 
Therefore, the first stage of the initialization procedure is to find the static 
solution p of the mean field problem at zero temperature with a small exter- 
nal constraint Ar 2 . This provides us with a well defined one-body potential 
U [p ] . Then, we have constructed a finite-temperature particle distribution 
as a hot gas of independent particles in this self-consistent potential, U [p ]. 
Then using the scaling 

r — > r' = ar, (3) 
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which conserves the entropy, we can generate various families of isoentropic 
initial conditions with different dilutions or compressions. 

2.4 Alternative approximate isothermal initialization 

If we do not take into account the self-consistent part of the mean-field po- 
tential, we see that energies scale approximately as 

e — ► e' ~ e/a 2 , (5) 

showing that the temperatures approximately depend on a as 

T — ► T' ~ T/a 2 , (6) 

Therefore we can also generate a family of approximately isothermal ini- 
tial densities by keeping T' constant, i.e. changing the temperature T prior to 
the scaling. Results obtained with this approximation are in all points com- 
parable to Constrained Hartree-Fock (CHF) simulations for a large range of 
initial conditions within the limitations of CHF discussed above. This gives 
good confidence in the scaling method and in the robustness of the conclu- 
sions we draw. Moreover, we would like to stress that, in fact, any initial 
conditions can be studied since the most important point is to use identi- 
cal protocols in quantum and classical simulations so that the comparison is 
meaningful. 



3 Quantum model 



Following fll3[ ) we consider spherically symmetric nuclei. The one-body wave 
functions can be separated into their radial part, angular part and spin- 
isospin part 

* a (?,s,t) = ^-Y lm (9,^ Xa>T (s,t) (7) 

in this expression f — (r,9,<p), x UT represents the spin and isospin part. In 
our special case, we consider that our system is spin and isospin saturated. 
Finally, a is an abbreviation for all the quantum number (n, I, m, a, r). 
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For the considered spherical system, the TDHF equation reduces to a 
radial equation for the single particle wave functions 

(8) 

In this equation V cou i om b is the Coulomb potential^ while U is the density- 
dependent mean-field potential 

U [p] (r, f) = ^oP(r, + ^V +<J (r, 0, (9) 
In this expression, p is the local density 

p(r,t) = 4H21 + l)^ l ^ )|2 , (10) 

n,l 

where n n \ are the time-independent occupation numbers which only depend 
upon n and I quantum numbers. The factor 4 is coming from the spin and 
isospin degeneracy 



3.1 Quantum initialization 

As discussed in the introduction, we have studied different types of initial 
conditions consisting in varying the dilution, the temperature and the entropy 
of the considered nuclei. 



3.1.1 Constrained Hartree-Fock 



In the constrained Hartree-Fock approach, the occupation numbers n n i are 
distributed according to the Fermi-Dirac distribution 



n n i{T) = 




where the energies e n i are the eigenvalue of the constrained Hamiltonian h + 
Xr 2 while the single particle wave functions are nothing but the eigenstates of 
the constrained Hamiltonian. The parameters A and \i are defined requiring 
a given dilution and a given number of particles. 

2 In this calculation the coulomb field is considered in an approximate way by giving 
the charge +1/2 to each nucleons. 
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3.1.2 Scaling initialization 



In the second method, we use the energy levels and the wave functions of 
the zero temperature mean field in presence of a small external field (Ar 2 , 
A = 0. 25 M eV fm~ 2 ) . This gives the energies e n i and the associated single 
particle wave functions R n i . Then, these states are occupied according to 
the Fermi-Dirac distribution ([H]). The temperature, T, is computed in order 
to get the desired entropy 

S = -4^(2/ + 1) (n nl \ogn nl + (1 - n nl ) log(l - n nl )) (12) 

n,l 

Finally, in order to get the appropriate dilution, the wave functions are 
rescaled by the transformation 

Klir') = ?0, (13) 

This initialization procedure can also be used in order to generate scaled 
density at approximately constant temperature by scaling the temperature 
T according to Eq. (IB). 



3.2 Numerical details 

In practice, the self-consistent solution of the static Hartree Fock (at fi- 
nite temperature and with a constraint Ar 2 ) has been found through an 
iterative application of the operator 1/ (e n i — (h + Af 2 ^ to the single par- 
ticle states In this procedure, at each step the one-body density is 
computed using the occupation numbers ( |TTj) where the energy is given by 
e n i = ($ni \ h + Af 2 |$„z). We start with the oscillator wave functions and 
after many iterations, we converge to the solution of our self-consistent field. 
We have controlled that this method gives same results as an iterative di- 
agonalization of h + Af 2 . From a practical point of view, we have done the 
simulation with up to 40 orbitals. We have chosen a step of Ar = 0.2/m in 
r-space. 

After having initialized the system, we let it evolve according to the 
TDHF equation written in r-space(|§|). During the expansion of the system, 
occupation numbers are kept fixed, implying the conservation of the entropy 
of the total system. Note that, in the following, we will sometime study only 
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part of the system; for such a sub-system the evolution will not be isoentropic 
anymore. 

We have performed simulations on a lattice of size 480/m, this size is 
big enough to avoid bouncing of particle against the wall. For the evolution, 
we have taken a step in time of 0.75 fm/c. Before presenting results, we 
will present in an analogous manner the method we have used for the semi- 
classical simulation. 



4 Semi-classical model 

The equivalent semi-classical equations are obtained by taking the first term 
of the h expansion of Wigner transform of the TDHF equation. This leads 
to the Vlasov equation 

9 -^ = W],f] (M) 

where {., .} are the usual Poisson brackets and f(r,p,t) and h[f] are respec- 
tively the one-body phase-space density and the associated self-consistent 
mean-field. One of the methods widely used for solving the Vlasov equation 



is the lattice Hamiltonian method developed in ref. |16 . 

In this method, the r-space is discretized into a three-dimensional cubic 
lattice. Each nucleon is represented by a number N test of "test particles". 
The density is calculated on each site of the lattice using the following algo- 
rithm: considering a site " i" , each particle contained in the neighboring cubes 
participates to the density at the point " i" with a given weight which depends 
on the distance between the particle and the considered site. Formally, this 
is obtained through the convolution 

p = p®G (15) 

where the unfold density is simply given by the density of point-like test- 
particles 

I ANtest 

p(x i ,y i ,z i ) = — — 5(xi - x n )5(yi - y n )5(zi - z n ) (16) 

latest n= i 

and where the weight function is given by a distribution of the density over 
2m neighboring sites. 
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G(x,y,z) = g(x) g(y) g(z) (17) 

, . (mAr — \a\) n , . . 

g (a) = , 7 ' - 9 mAr- a) 18 

In this expression, corresponds to the Heavyside function, so that (mAr) 
corresponds to the range of the smoothing used to compute the density. 
Finally we get the following density at each site i of the considered lattice 

-y ANtest 

p(x h Ui, Zi) = — — 9(xi - x n )g(yi - y n )g(Zi - z n ) (19) 

^test n=1 

This density is used to compute the energy functional E [jo] (see eq. (0)) at 
each point i of the lattice, in order to define the Hamiltonian 7i of the system 

2 

n(f n , Pn) = E 1^ + E Ar " E [PI (^) + Ecoulomb (20) 
n ^ m i 

where, following our conventions, the sum over n runs over all the test- 
particles while the sum over i corresponds to the sum over lattice sites. The 
evolution of the system consists in the classical Hamilton equations for each 
test particle n which have been solved using the standard leap-frog algorithm 

f(t + At) = f(t - At) + 2At — , (21) 

m 

p{t + 2At) = p{t) + 2At F(t + At), (22) 
where the force is calculated according to 

- dH(f n ,Pn) 



df r , 



(23) 



Using the expression (|20"D for the Hamiltonian 7i, it is easy to demonstrate 
that the force F derives from a potential U which is related to the mean field 
potential U [p] through the relation 

U (x, y, z) = -^—^ E U ( x *> Vi> z i) ( 24 ) 

where the sum runs over the 2m neighboring cells of the position (x, y, z). 
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This method has demonstrated its accuracy and is now widely used in 



theoretical description of Heavy- Ion collision [16]. We have thus used it and 



compared results with quantum simulation. It should be noticed that other 
methods are reported in the literature using either slightly different lattice 
discretization of the density or gaussian smoothing in order to fold the point- 
like test-particle density into a smooth nucleon density. Therefore, all the 
semi-classical methods contain a smoothing procedure in order to generate 
the nucleon density used in the energy functional. The presented results do 
not depend on the particular numerical implementation of this smoothing. 

In practice, we have adopted the value N test = 200 which corresponds to 
a rather large number in such type of simulations. For the Vlasov evolution, 
we have taken the three dimensional code TWINGO developed in ref. [TBI ITS 



The simulation are done with a range for the convoluting function g of 2 fin, 
the lattice grid has a step of lfm and the step in time is the same as in the 
quantum simulation: 0.75/m/c. The total r-space is a box of size 40 x 40 x 40, 
when a particle reachs the boundary, it is not evolved anymore. 

Before describing the results, we will get more insight into the simulation 
and in particular into the initialization procedure. 



4.1 Initialization and stable states of the Vlasov dy- 
namics. 

4.1.1 Constrained Thomas- Fermi solution 

The first point is to define stable spherical nucleus at temperature T under 
the constraint Ar 2 . For this purpose, we have only to compute the spherical 
one-body phase space density f(r,p) of the Thomas- Fermi type 

f{r,p,T) = I exp I I +11 (25) 

which indeed maximizes the entropy 

S = aJ.-.J (/log/ + (1 - /) log(l - /)) (26) 

under the various constraints. In expression (f25j), /i is the chemical potential 
and U(r) is the potential consistent with the lattice Hamiltonian method (see 
Eq. (p4|)) which is obtained using the folded density p. Note that, in order to 



11 



speed up the initialization process, we have used the spherical symmetry: we 
have replaced Eqs. ([19]) and (|24f) by convolutions with a spherical function 



p(r) = / p( r ') g( r ' ~ r)Airr' 2 dr' (27) 

Jr=0 

_ rr'=+oo 

U(r) = / U[p]{r') g(r' -r)A-Kr' 2 dr' (28) 

Jr=Q 

where the function g is a function analogous to the expression fll8|) 

^) = ^r^e(r -r) ( 29 ) 

but with a range r adapted to provide the same averaged squared-radius as 
G. 

The numerical procedure used to build a constrained Thomas-Fermi dis- 
tribution in phase space can be separated in two steps: 

• the first one consists in an iterative definition of f(r,p, T): starting with 
a potential (at the first iteration the potential is assumed to be a Wood- 
Saxon), we calculate the phase space density according to expression 
(p5|) using a chemical potential p which conserves the particle number. 
Then, we obtain the density p(r) by integrating f(r,p, T) over p and we 
calculate the folded density p{r) and the associated potential U (r) using 
Fast Fourier transform techniques. Finally, we iterate the operation 
until a full convergence is reached. This provides us with a continuous 
self-consistent phase-space density f(r,p,T). 

• The second step consists in distributing test particles according to the 
phase space distribution, f(r,p,T) using a standard Metropolis algo- 
rithm [ 17| . 



The stability over time of a nucleus initialized in its ground state demon- 
strates the accuracy of this method in finding stable solutions in lattice 
Hamiltonian framework.^ 

3 Indeed, we have not observed the emission of a single test-particle during the typical 
simulation time used in the following, (more than 300/m/c) 
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4.1.2 Scaling procedure 

As already discussed in the quantum case, we have also studied a second 
class of initial conditions namely those obtained through a scaling. In order 
to follow the approximation made in the quantum simulation, the mean- 
field potential, included in the definition of the finite temperature phase 
space density (eq. Q2"5D), has been approximated by the zero-temperature one 
(including a small constraining external potential Ar 2 , A = 0. 25 M eV fm~ 2 ). 
For each of these initial conditions, the temperature is chosen in order to get 
the expected entropy (|26|).Then, the test-particle distribution is generated 
using the Metropolis algorithm as discussed above. Finally the momenta 
and the positions of the test particles are scaled according to eq. (|3|). Since 
this scaling preserves / and the phase space volume d s rd 3 p/h 3 , it does not 
change the entropy (p6|). We are thus generating a family of isoentropic 
initial conditions. 

This scaling procedure can be also used to construct an ensemble of initial 
conditions which can be considered as isothermal by changing the tempera- 
ture T in order to compensate the effects of the dilution (or the compression) 
on the temperature prior to the scaling according to Eq. (|). 

4.2 Discussion on smoothing in semi-classical methods 

In this section, we will discuss the possible bias of the comparison between 
quantum and semi-classical simulations due to the test-particle method. 
Indeed, in any semi-classical approach based on a test-particle method, a 
smoothing of the phase space density should be included in order to com- 
pute the mean-field potential. As already mentioned in PS[, this smoothing 



acts as an effective range for the forces. Indeed, considering the density of 
particles p, the actual numerical smoothed density p is given by the folding 
product with a smooth function g 

p = p®g (30) 



The function g used in our simulations is given by equation pq ) but g could 
take other expression. Expending the Fourier transform of g at the second 
order in k this folding product can be approximated by 

p = p-f3Ap (31) 
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where (3 is the range of the function g. Then, the time- dependent mean-field 
becomes 



U [p] = f p + -^^P +a ~ P A P (j*o + - ^ L hff\ (32) 

where one can recognize a typical finite range term in the last part of Eq. 



(|32|). However, it should be noticed that the above expression corresponds 
to a very specific surface contribution. 

As a direct consequence, due to this additional surface energy, we will 
see that energies calculated in semi-classical approximations will be greater 
that those calculated in quantum simulations. This might seems surprising 
since in absence of quantum correction terms in the semi-classical potential, 
one would expect the surface energy present in the semi-classical picture to 
be smaller than the quantum equivalent (or even zero). In fact, this surface 
energy is only generated by the numerical method. Indeed, in the quantum 
case, since wave functions are not localized, no smoothing is needed and the 
algorithm does not introduce spurious additional surface term. However, it 
is well known that even with a zero range force with no momentum depen- 
dence (to simulate a finite-range), a quantum calculation will exhibit non 
trivial surface properties due to the role of the kinetic term and because of 
the self consistence effects on the wave functions. Conversely, a semiclassi- 
cal calculation with a zero range force should present sharp edges (at zero 
temperature). Then, surface properties, in semi-classical calculation, are in- 
fluenced by the smoothing procedure. 

This difference in the surface properties can be seen as a major drawback 
of semi-classical approaches. Of course one might think to take advantage 
of the presence of this numerical finite range in semi-classical pictures in 
order to mimics some quantum effects or some surface terms in the consid- 
ered interaction. In the next chapter, we will more carefully investigate this 
possibility, however, we can already mention some difficulties: 

• It is never easy to rely on the numerics in order to mimic physical 
effects. 

• One might try to select the numerical coarse-graining, but normally, it 
is fixed by some other arguments such as numerical accuracy (see Eq. 
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since the surface effects are coming from the folding procedure in the 
mean-field potential their importance is fixed by the mean field strength 
(see Eq. (|32|)). Since this strength is density dependent, the surface 
effect will have the same dependence. This might be not realistic since 
they are supposed to mimic surface effects of various origin and so with 
various density dependence: interaction range, quantum kinetic term, 
finite range of wave functions and self-consistent effects. This point is 



partly discussed in ref. [23 



In any realistic cases, where surface term exists in the interaction, one 
should be careful in order to not double count these finite size effects. 



We observe that, a semi-classical simulation, which is supposed to be 
based on a given force ([[]), includes a spurious additional interaction ([!2|). As 
we will see, this will contribute to the observed differences with the quantum 
simulation associated to ([!]). However, realistic interactions usually include a 
finite range term, so that it might be interesting to investigate the possibility 
to mimics quantum dynamics using a finite range interaction with a semi- 
classical approach eventually taking into account the numerical surface term 
in the potential. 



4.3 Discussion of surface effects. 

In this section we consider a mean-field potential taking into account a finite 
range effect in the TDHF mean-field^ 

U( P ) = 2 f„ P + h P° +1 + c VV (33) 

Usually, in the TDHF calculation, the Laplacian terms is simulated by a 
finite range potential of the Yukawa type in order to avoid well-known numer- 
ical instabilities |24]]. Therefore, the mean-field potential has been computed 
using the relation 

U(p) = \kP+ h P° +1 + V oP ®Y. (34) 

4 In the skyrme force, the parameter c is related to the parameters t\ and £2 Hi PJJ- 
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where Y is a Yukawa folding function 

nr) = ^£|ZMM. (35) 

and where t' is function of a, Vo and the to used in the classical calculation. 
Indeed, expanding the Fourier transform of the folding product up to the 
second order in k we can derive the following expression for the potential 



U(p) = 3 -(t + p + ^±2) t3 + 2W V V (36) 



For numerical applications, we have chosen the same values of to and £3 as 
before in the semi-classical simulation. For the quantum case, we have chosen 
for a and Vo respectively 0.37 fm and —200 MeV (which gives a t value of 
830.3 MeV fm 3 ) in order to have almost the same energy for the ground state 
in quantum and in semi-classical case at zero temperature. In the following, 
this particular simulation will be referred as " Quantum+ surface" . 



5 Results 

In this section we will focus on a detailed comparison of the two considered 
approaches, the semi-classical Vlasov dynamics and the full quantum TDHF 
treatment (with and without additional surface terms). We will take a hot 
and compressed or diluted Double magic 40 Ca nucleus as a test case. 

5.1 Initialization of nuclei under extreme conditions 

In this section, we will illustrate different initial conditions of compression 
and temperature. We will first discuss density profiles obtained after the 
initialization step and then, present the various excited nuclei we are consid- 
ering in the evolution. 

5.1.1 Density profiles 

In figure (|1]) , we compare quantum and semi-classical constrained mean-field 
solutions obtained at various temperature (T = 0, 5, 10 and 15 MeV) with a 
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small constrained field Ar 2 with A = 0.25 MeV/fm 2 . As expected, the semi- 
classical simulation, which, in addition to the smoothing procedure, does 
not include shell effects, presents smoother initial density profiles. Directly 
connected to the smoothing, the quantal density profiles always present a 
smoother surface than the semi-classical one. In addition, the Root Mean 
Square radii (RMS) and the central densities are reported in table [T| We 
observe that these quantities are slightly different in semi-classical case com- 
pared to both quantum simulations. In particular, the temperature depen- 
dence of such global observables indicates a different response of the nuclei 
to the heating process within the two approaches. 




0246024602468 



r(fm/c) r(fm/c) r(fm/c) 



Figure 1: Initial density profiles of a 40 C a nucleus at various temperatures 
for a constraining field Ar 2 with A = 0.25 MeV/fm 2 : a)Quantum approach; 
b) Quantum+surface c) Semi-classical. In each figures are displayed temper- 
ature T = 0MeV(solid line), 5MeV (dashed line), lOMeV (long-dashed line), 
15MeV (short-dashed line). 
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Classical 


RMS 


P 
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P 


T = OMeV 


3.14 


0.172 


3.31 


0.153 


3.09 


0.196 


T = 5MeV 


3.30 


0.148 


3.53 


0.132 


3.36 


0.179 


T = lOMeV 


4.30 


0.130 


4.65 


0.109 


4.93 


0.144 


T = 15MeV 


5.55 


0.098 


5.94 


0.082 


8.03 


0.080 



Table 1: Root Mean Square radius (RMS) and central densities of a 40 Ca 
nucleus initialized with a constraining field Ar 2 with A = 0.25 MeV/fm 2 at 
various temperatures 

In the scaling procedure, for isoentropic initial conditions, density profiles 
of compressed or dilated nuclei will be obtained by a scaling of those displayed 
in Fig.0. 

5.1.2 Isoentropic initial conditions 

Important information about a system is provided by the total energy E as a 
function of p for various entropies. When considering infinite nuclear matter 
this is the usual Equation Of state (EOS), E = E(p, S) as shown in fig. 2 for 
the considered interaction. This EOS was extensively used in order to under- 
stand the properties and the dynamics of nuclear matter. For each entropy, 
the minimum of the energy corresponds to a stationary isoentropic solution. 
Furthermore, the second derivative of the energy with respect to p around 
the stationary point is directly related to the frequency of the monopole 
oscillation around this point. The isoentropic spinodal region is defined as 
points where this curvature is negative or equal to zero. This region is of 
particular interest since one of the possible scenarii of multifragmentation is 
that a highly excited nucleus could enter this region and become unstable 
against the partition of the matter in fragments. 

Let us built an equivalent of the EOS for a finite system. In order to do 
so, we can plot the relation between the central density p and the energy per 
particles E/A of the system for the various initial entropies considered. In 
Fig.^|, we report four different curves which correspond to entropies S/ks = 
0, 1.10, 2.35 and 3.28. These curves are obtained by applying the general 
scaling procedure to each models: occupations numbers are kept fixed while 
different compressions/dilations are considered. 
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Figure 2: Binding energy of the infinite nuclear matter as a function of its 
density, the curves correspond respectively to S/k B = 0, 1.1, 2.35 and 3.28. 
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Figure 3: Binding energy of a hot and compressed/diluted Ca nucleus as 
a function of its central density, the curves corresponds respectively from 
bottom to top to S/k B — (x), 1.1 (o), 2.35 (□) and 3.28 (o) on each figures. 
Left: Quantum (without surface). Middle: Quantum+Surface. Right: Semi- 
classical. 



20 



Looking first at the minimum of the curve S/ks = 0, we see that the 
semi-classical stationary point has an energy greater than that of the quan- 
tum system (without surface). This is a direct consequence of the smoothing 
procedure. The additional surface term gives an extra energy of the order 
of 3 MeV. Note that, in the second quantal calculation (middle part), the 
surface term was adjusted in order to have the same ground-state energy as 
in the semi-classical case. If we now focus on the curvature of isoentropic 
curves, we see that they appear smaller in semi-classical than in both quantal 
calculations. Note that, if we suppose an evolution at contant entropy and 
energy, these curvatures are direcly related to the response to a monopolar 
excitation. In such a context, according to figure the semi-classical model 
should have a larger amplitude than quantal one. However, in dynamical 
simulations, physical systems have neither constant energy nor constant en- 
tropy in particular due to particle evaporation. As we will show later, the 
difference in evaporation processes in addition to the damping due to the 
use of test-particle method will strongly reduce the amplitude of monopolar 
vibrations in semi-classical simulations making it smaller than the quantum 
analogue. 

Finally, we want to point out that temperature creates more dilution and 
disorder in semi-classical calculation than in quantum one as we already no- 
ticed in Fig.|I|. We can quantify this more easily considering for instance 
the dependence of the entropy with temperature at a fixed compression fac- 
tor. In Fig. ^, we report such a dependence for the three models when no 
compression is applied. We see in particular, that higher temperatures are 
necessary in quantal cases in order to create the same entropy than in a 
semi-classical simulation. This is a great difference between quantum and 
classical treatment of nuclei. 

5.2 Monopole oscillations and thermal expansion 
5.2.1 Examples of expansion 

Let us now consider the time evolution of the density profiles presented in 
figure [I] for the constrained method. In figure |^, we present the evolution 
of a hot system for different temperature T = 5, 10, 15 MeV. From a 
qualitative point of view, the semi-classical and the TDHF evolutions are 
in reasonable agreement except at high temperature, for which, the semi- 
classical nucleus does not resist to the heat and to vaporizes easily than 
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Figure 4: Entropy as a function of the temperature for a 40 Ca nuclei obtained 
when a small external field Ar 2 with A = 0.25 MeV/fm 2 is applied. Solid line: 
Quantum (without surface), dashed line: Quantum+ surface, long-dashed 
line: Semi-classical. 
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its quantum analogue. In fact, not only the vaporization process appears 
different and seems to start at a lower value of the initial temperature but 
also the whole particle evaporation is faster. 

Before investigating in more details this observation and before discussing 
the possible reasons of this behavior let us first study quantitatively the var- 
ious observed differences between the semi-classical and the quantum ap- 
proaches. 



5.2.2 Illustration of the monopole oscillation 

Let us first study the collective monopole motion associated with a com- 
pressed or diluted system. 

Figure ^ shows the central density as a function of time, predicted by 
the quantum evolution for a heated and compressed Ca nucleus initialized 
at T = 5 MeV and with an initial density p ini = 0.23 fm -3 . In all cases, the 
central density presents an oscillating pattern due to the monopole vibration 
of the nucleus. In addition, in Quantum simulation beatings of different 
modes which are characteristic of the Landau spreading of the breathing 



mode onto various components [21, 22], are observed even at high temperature 



(see Fig. [5] (Top)). Note also, that the central density in semi-classical 
simulations , exhibits an important noise added on top to the monopole 
vibration. These fluctuations are due to the finite number of test-particles. 
In order to quantitatively study this vibration, for different initial conditions, 
we have considered two observables as a function of the initial entropy and 
density: 



• the value of the minimum of the central density reached during the 
evolution (noted p min which provides a measure of the amplitude of the 
monopole vibration; 

• the time needed to reach this minimum (noted i m i n ) which, in case of a 
pure oscillatory motion, is directly related to half of the period of the 
vibration. 



In the following, we will present results where events are grouped accord- 
ing to their entropies. 
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Figure 5: Density profiles for various initial temperatures - for T = 5, 10, 15 
MeV (from the top to the bottom) - and for: Left part: Quantum (Without 
surface). Middle Part: Quantum+Surface. Right part: Semi-classical. In 
each figure various curves corresponds to different times: Solid line: t — 
fm/c. Dashed line: t = 30 fm/c. Long-dashed line: t = 300 fm/c. 
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Figure 6: Illustration of oscillation of the central density in time for different 
models. In all models, initial conditions correspond to a 40 Ca nuclei at T = 5 
MeV and an initial density p ini = 0.23 fm~ 3 . Top: Quantum (Without 
surface). Middle: Quantum+ surface. Bottom: Semi-classical. 
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Figure 7: The central density at the maximum of the dilution (turning 
point) as a function of the initial density for a 40 Ca nucleus at various ini- 
tial temperatures (from top to bottom points corresponds respectively to 
S/k B = (x), 1.1 (o), 2.35 (□) and 3.28 (o)). Left: Quantum (without 
surface). Middle: Quantum+ Surface. Right: Semi-classical simulation. 
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5.2.3 Expansion of isoentropic initial conditions 

The value of the first minimum of the central density oscillations (p m j n ) is 
displayed in figure (0) as a function of the value of the initial density (p init ). 
For a given entropy, this curve presents two branches. The first one (low 
initial density region) corresponds to the case of a minimum density equal to 
the initial one, i.e. the nucleus is already diluted compared to the stationary 
initial condition. Therefore this first branch corresponds to p min = p init . The 
second branch gives the amplitude of the collective motion, i.e. the maximum 
dilution, as a function of the initial compression. The point where these two 
branches met is the stationary point for which no oscillations are observed. 
Note also that at finite temperature even if the system is stationary, it cools 
down by particle evaporation. 

Comparing now the classical and quantal simulations, one can see several 
differences. First of all, it is clear that the effect of entropy (or temperature) 
is weaker in semi-classical simulations when the excitation is not too high. 
In particular the stationary points do not exhibit the same entropy (or tem- 
perature) dependence in the quantum case or the classical case (see. Fig. 

In a similar manner, the initial densities for which the central density will 
eventually reach (i.e. either the nucleus vaporizes or leads to an hollow sys- 
tem) have different temperature dependence (see Fig.^J) . At low temperatures, 
semi-classical approaches are associated with a fast emission of particles. This 
faster evaporation as compared with the quantum dynamics induces a faster 
damping of the collective motion. Moreover, it is known that the fluctuations 
induced by test particle methods as well as the possible turbulent dynamics of 
classical fluids generate an additional damping of collective motions in semi- 



classical approaches pq| . All these different dampings reduce the amplitude 
and the frequency of vibrations. In such a way, a more compressed nucleus 
can survive more easily in semi-classical simulations because it can dissipate 
energy more efficiently than its quantum analogue and has a longer time to 
cool down. At high temperature, all particles are immediately evaporated in 
the semi-classical approaches leading to the vaporization of the nucleus. 

If we now look at the slopes of p min versus p ini in the right part of figure [7|, 
we can see that the slopes appear smaller in semi-classical case than in both 
quantum cases. This behavior indicates a larger damping of the collective 
motion together with a softer response to an external compression of the 
nuclei in the semi-classical treatment, compared with the quantum ones. This 
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Figure 8: Stationary points obtained from the (p min , Pmit) diagram. The 
stationary density value is plotted as a function of the entropy. Solid-line: 
Quantum (without surface). Dashed-line: Quantum +Surface. Long dashed 
line: Semi-classical. 
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Figure 9: Initial density for which the central density reach during the 
evolution (from the (p min , Pi n it) diagram). The initial density value is plotted 
as a function of the entropy. Solid-line: Quantum (without surface). Dashed- 
line: Quantum+ Surface. Long dashed line: Semi-classical. 
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apparent softness could appear surprising a priori since we are considering 
models with same EOS for the infinite nuclear matter. However, we have 
seen in the previous section that the EOS for finite systems are different in 
the various simulations (see Fig. |3|). 

In conclusion, a higher compression is needed in semi-classical treatments 
as compared with quantum approaches in order to reach the same minimum 
density. This corresponds to a smaller amplitude and a slower collective 
motion. Part of this difference could be explained by the smoothing proce- 
dure necessary in semi-classical method to generate the mean-field. How- 
ever, adding a surface term in quantum mechanics reduce only a little this 
difference, indicating that physical aspect such as, for instance, dampings, 
evaporation, time-scale... are different. In the following, we will quantify 
differences in time-scale and in the evaporation process. 

5.2.4 Monopole vibration period 

In the previous section, we have mentioned that semi-classical methods lead 
to much softer response of the nucleus. As a direct consequence, we could 
expect a slower expansion dynamics in the semi-classical treatment. In Fig. 
PH| , we have plotted the time (noted t min ) necessary to reach the lowest central 
density p min for isoentropic initial conditions. At high initial compression, 
this time goes to infinity since the system is vaporized. For initial dilatation 
below the stationary point, this time is simply equal to zero (or nearly zero) 
since the system immediately contracts. Inbetween, we have finite, non- 
zero t min which could be interpreted to be one-half of the period of monopole 
vibration. In the two quantum cases, at temperature or entropy equal to zero, 
the monopole vibration period is close to 80 fm/c and appears to be almost 
independent of the initial perturbation (which is similar to the harmonic 
picture). In the classical case, the period of monopole vibration is of the order 
of 120 fm/c. This indicates again a softer response of the classical system 
accompanied with a longer time of expansion. If we now look to higher 
temperature (or entropies), we see that in quantum cases, the monopole 
motion period increases with temperature (or entropy) and that for a given 
temperature or entropy, this time increases with the initial compression. The 
situation is not as clear in the semi-classical case mainly due to the presence 
of large fluctuations in the central density which makes the extraction of the 
corresponding times difficult. 

Therefore, as far as the collective motion is concerned, we can conclude 
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Figure 10: Time t min to reach the minimal density as a function of the initial 
density p init at various initial entropies S/ks = (x), 1.1 (o), 2.35 (□) 
and 3.28 (o). Left: Quantum (without surface). Middle: Quantum+Surface. 
Right: Semi-classical simulation. Note that, near the stationnary point, the 
time t min could not be extracted in semi-classical simulation due to the strong 
fluctuations induced by test-particle method (see figure 6). These points are 
thus not represented in the figure. 
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that, compared to the semi-classical predictions, the quantum simulations 
present a faster expansion towards lower densities. These differences appear 
to be due to the differences in the collective potential but also to the dif- 
ferences in the damping and cooling processes which appear faster in the 
semi-classical case. Then, in order to conclude this comparison let us discuss 
in more details the evaporation process. 

5.2.5 Evaporation dynamics and time scales 

In order to study the evaporation process we have computed the particle 
flows outside of a sphere of radius 15 fm. These flow are shown on figure 

According to this figure it is clear that the semi-classical approximation 
predicts a quicker evaporation of more particles than the quantum one. We 
see that these differences in the evaporation process are reduced by the in- 
troduction in the quantum simulation of a surface interaction (see right part 
of Fig.|ll]). This can be clearly quantified by the asymptotic number of evap- 
orated particles (after 300/m/c) which is presented in Fig. ^ as a function 
of the temperature/entropies for different simulations. We can see that the 
temperature acts completely differently in quantum and semi-classical cases. 
In particular for high temperature, the semi-classical model is poor approx- 
imation of the quantum mean-field theory. The heated nucleus appears less 
bound in semi-classical than in quantum cases. However, we note that the 
introduction of a surface term reduces quantitatively the gap between the two 
treatments. Nevertheless, the time scales of the evaporation are different. 

This important difference can be traced back to the fact that in a quan- 
tum approach the nucleons can be partly reflected by the mean-field potential 
wells while in a classical case as soon as a particle as enough energy it will be 
immediately emitted. This large difference between semi-classical and quan- 
tum approaches is one of the major drawback of semi-classical approximation. 
Indeed, it affects not only the evaporation itself but in fact the whole dynam- 
ics changing time-scales, cooling processes, influence of the thermal pressure 
and also the one-body dissipation. 
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Figure 11: Particle Flow out of a sphere of radius 15 fm for a typical Ca 
nucleus at various temperatures T = 5, 10 and 15 MeV (from top to bottom). 
Left: Quantum without surface (dashed line). Right: Quantum+Surface 
(dashed line). On each figures is superimposed the Semi-classical case (solid 
line). 
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Figure 12: Asymptotic number of evaporated particles (after 300/m/c) as a 
function of the temperature for the different considered models. Solid line: 
Quantum without surface. Dashed line: Quantum+Surface. Long-dashed 
line: Semi-classical. 
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6 Discussion of the expansion and spinodal 
instabilities 



Semi-classical simulations (including fluctuations) of the spinodal decompo- 
sition are able to reproduce globally the experimental fragment partitions [fT9|, 
However, these simulations seem to underestimate the kinetic energies of 



the fragments. This is a general observation made on BUU simulations that 
the expansion velocity at the entrance of the spinodal region is barely enough 
compared with the experimental data leaving small room for additional slow- 
ing down during the fragment formation. It is thus important to discuss the 
differences observed between semi-classical and quantum dynamics when the 
nuclei reach low density regions. 

At zero temperature, we have seen that the infinite system becomes un- 
stable when the density is less than O.l/m -3 . 

From figure [13|, it is clear that higher densities (up to a factor 2) are 
needed in semi-classical calculations to reach a density below 0.1 fm~ 3 . If 
we look now at figure O, we see that the time needed to reach the low 
densities is also longer in semi-classical simulation. This should be related 
to the observed stronger damping of collective motion in the Vlasov case. 
The induced delay of a factor 1.5 to 2 is associated with a reduction of 
the collective expansion velocity by the same amount. This may explain 
why semi-classical simulations taking into account the fragment formation 
dynamics seem to underestimate the experimental kinetic energies. From 
this study, it is now clear that quantal simulations are important when one 
wants to discuss quantitatively collective expansions and fragment formation 
that could occur during Heavy Ion reactions. 



7 Concluding remarks 

This work presents a comparison between quantum and semi-classical mean- 
field dynamics of hot and compressed nuclei. 

We have first pointed out some properties of the semi-classical approaches 
using test particle sampling of the phase space. One particularity of this 
method is the smoothing or the coarse-graining of the phase space in order 
to define the density. This smoothing adds a finite range term to the initial 
potential (expression (|313)). This numerical surface term affects static and 
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Figure 13: Initial density that gives a minimal density equal to 0.1 fm -3 as a 
function of the entropy. Solid-line: Quantum (without surface). Dashed-line: 
Quantum+ Surface. Long dashed line: Semi-classical. 
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Figure 14: Time necessary to reach a minimal density equal to 0.1 fm~ 3 as a 
function of the entropy (for the particular initial density that has the minimal 
density at 0.1 fm~ 3 . Solid-line: Quantum (without surface). Dashed-line: 
Quantum+Surface. Long dashed line: Semi-classical. 
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dynamical properties of nuclei. We have disentangled effects due to this extra 
surface term from intrinsic differences due to the absence of quantum features 
in semi-classical simulation by comparing them with quantum calculations 
in which a short range potential was added in order to mimic the numerical 
one present in semi-classical approaches. 

From these three different types of simulations (semi-classical, quantum 
with and without finite range interaction) we have drawn the conclusions 
listed bellow. 



While the three approaches have identical static properties of infinite 
nuclear matter, the energies of finite size nuclei at different dilutions 
and excitations are different. In particular, the compressibility modulus 
(or equivalently the curvature of the EOS) of the finite systems depends 
upon the considered approach. This directly affects the collective dy- 
namics in nuclei. In particular, we have shown that the frequencies of 
the breathing mode, the dependence in temperature of the amplitude of 
the collective modes and the position of stationary points are affected 
by the presence of a numerical gradient of density in the potential. 
However, the presence of this numerical surface term do not explain all 
the observed differences and part of the discrepancy is also due to the 
semi-classical approximation. 

From the dynamical evolutions, we can also see that the damping of 
the collective mode is greater in semi-classical method than in quantum 
approaches despite the absence of Landau spreading of the breathing 
mode []2T| in semi-classical approaches. This strong damping observed 



in semi-classical simulations reflects the differences in the evaporation 
process as well as additional dissipation in the classical motion of par- 
ticles. In the semi-classical case, highly excited particles are emitted 
early in the evolution inducing a large energy loss during the initial 
stage of the dynamics. This implies a fast attenuation of the collec- 
tive mode. Conversely, in the quantum case, the reflections of the 
single-particle wave functions on the potential well reduce the particle 
evaporation rate and increase the pressure on the nucleus boundary. 
This difference on time scale is smaller but remains important when a 
surface interaction is introduced in the quantum-treatment. The evap- 
oration is faster and the total number of emitted particles is larger in 
semi-classical treatments compared with quantum approaches. 
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• The differences in the evaporation process not only affect the damping 
of the monopole vibration but they also change the size of the remain- 
ing residue after desexcitation. At low temperature but extreme con- 
ditions of compression, the fast evaporation observed in semi-classical 
approaches gives a strong loss of energy that enable the system to sur- 
vive. On the contrary, at high energy, in the semi-classical case, the 
potential is not able to retain the particles inside the nucleus: the sys- 
tem is directly vaporized. In the quantum case, at low temperature, 
the collective motion is less damped so that the expansion leads more 
easily to a complete dilution of the system. Conversely, at high temper- 
ature, the evaporation is reduced by the barrier reflections so that a non 
compressed nucleus may survive the evaporation stage. In summary, 
in quantum approaches the mean-field is more robust. 

• Finally, we have discussed that in a semi-classical treatment the expan- 
sion is slower and that stronger compression are needed in order to reach 
low density regions. This induces an underestimation of the expansion 
velocity by a factor 1.5 to 2. This may fill the gap between theoreti- 
cal prediction of fragment kinetic energies and experimental data since 
up to now stochastic mean-field description of multifragmentation have 
been always performed using a semi-classical approximation. 
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